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ABSTRACT 

In this paper semi-active control of the suspension of an all-terrain vehicle is 
considered. A seven degree of freedom suspension model is presented first. A 
fuzzy approach for controller synthesis is then proposed. Expert knowledge is 
stored in the form of IF-THEN rules. The Takagi-Sugeno inference system is 
employed, with triangle membership functions. The fuzzy system output is the 
damper coefficient. In contrary to many other control algorithms, the presented 
fuzzy algorithm does not require inverse modelling of the magnetorheological 
damper. Instead, some scaling parameters are set. They can be chosen 
experimentally, or a bio-inspired strategy can be applied. The fuzzy control is then 
compared with the Skyhook control in simulations, in terms of road holding and 
driving comfort indicators. Obtained results are similar. However, lack of necessity 
to use an MR inverse model allows the fuzzy system to provide successful 
performance in case of different operating conditions, what is an important benefit. 
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I. INTRODUCTION 

Suspension of a vehicle, in addition to its mechanical role, when appropriately 
designed can successfully reduce vibration of the vehicle. Depending on applied 
vibration reduction method, three types of suspension can be distinguished: passive, 
semi-active, and active. Passive suspension does not use any extra force, and has 
invariant stiffness and damping. Its application results in significant vibration 
reduction. However, passive suspension parameters must provide a fixed-in-time 
trade-off between safety and driving comfort despite the driving conditions. In 
semi-active suspension, stiffness of that suspension can be changed during driving. 
Unlike passive suspension, semi-active suspension can be adapted to different 
conditions online. It can result in improving both safety and driving comfort. It 
involves a dramatically smaller energy amount than active suspension control, 
which requires actuators to generate additional force with a benefit of highest 
potential to reduce vibrations. 

Recently, growing interest in the use of semi-active suspension is observed, 
e.g. in [1], Such trend is caused by popularity of some materials such as 
magnetorheological (MR) dampers and their successful application to vehicle 
suspension control [2]. 

Skyhook is the most popular semi-active algorithm. Its idea is well explained in 
the literature, e.g. in [1]. There are also other semi-active algorithms, including: 
Groundhook [3], LQR [4], H°° Control [1, 4], and the group of smart algorithms 
including fuzzy control [5], neural networks [6], and genetic algorithms [7]. Other 
non-linear control algorithms, including adaptive properties, can also be adopted 
from the neighbouring active control field [8]. 
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This paper focuses on developing a fuzzy system to control semi-active 
suspension including a non-linear magnetorheological damper. The suspension is 
described in the form of a seven-degree-of-freedom model with parameters 
presented in [10]. Simulation analysis is performed based on data originating 
from a real all-terrain vehicle available for experiments. Results obtained with 
the fuzzy control algorithm are compared with those obtained with the Skyhook 
algorithm [11]. 

2. MODELLING 

There are many suspension and full-vehicle models in the literature. A significant 
amount of papers use 2 degree of freedom (2 DoF) quarter models [1, 4, 9], 
Recently, a 7 DoF model, also referred to as the ‘full-car model’, has become 
popular [3, 5, 7, 12, 15]. There are also publications that use different car models 
[10]. An 11 DoF model presented by Zareh, [10], includes additionally presence of 
four passengers. 

The chosen 7 DoF model reflects interdependencies between vibration 
transmission of vehicle’s quarters. Moreover, it provides the data about the roll and 
pitch angles, which are used to define the driving comfort. More analysis of vehicle 
suspension models can be found in [16], [17], [18], [19]. 

2.1. A seven degree of freedom passive suspension model 

The suspension model used in this paper is presented in Fig. 1. For a clear 
presentation of the model, Fig. 1 includes a full description for the unsprung masses 
m 2 , m 3 , and related displacements. 

For all simulations some assumptions are made, similar to those presented in [3], 
The road is assumed to be even, except for derbies. The whole construction is 
perfectly rigid. Parameters such as suspension damping, stiffness, and the tires 
stiffness are constant during the experiment. 

For the suspension model shown in Fig. 1, the dynamic equation can be 
expressed as: 

M ■ z M = ~[F U + F 2s + F 3s + F 4s + F w + F 2D + F 3D + F 4D ], (1) 

where: 

M -vehicle body mass, 

z M -acceleration of the centre of vehicle mass, 

F\ 4s , F j 4D -forces generated by springs, and dampers, respectively. 

The forces expressed in (1) can be defined as follows: 

F is = k iS (z is -z iu ), (2) 



Figure I. A 7 DoF suspension model. 
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(3) 


where: 

i -wheel index, 

k iS , c iD -spring stiffness coefficient, and suspension damper coefficient, 
respectively, 

z js , z iu , -sprunged and unsprunged masses displacement, respectively, 

z js , z iu -sprunged and unsprunged masses velocity, respectively. 

The dynamic equation for the unsprung mass acceleration is given as: 

m Au = F ts + f id ~ Kt (z iu ~ Zi ) (4) 

where: 

m i -unsprung mass, 

k iT -tire stiffness coefficient, 

Z t -road profile excitation. 

Torque equations for the roll cp and pitch 6 are defined as: 

Jfr<P = r [ F 2s + F 3s + F 2D + F 3o] COS <4) - l [ F ls + F 4s + F 1D + 4d] c0s (<?4 (5) 


J,j = f[F u + F 2s + F w + F 2d ]cos(0) - b[F 3s + F As + F 3D + F 4d ]cos(0), (6) 


where: 

J PS , Jpg -moment of inertia of the front-rear-axis, and the left-right 
-axis, respectively, 

f b, l, r -length of the chassis from COG to the front, rear, left, and 
right corners, respectively. 

2.2. The Bingham MR damper model 

In this paper the RD-1005-3 MR damper, with 0-5V input limit is used. For its 
modelling the Bingham model is chosen. Its structured is relatively simple and MR 
damper properties are sufficiently well represented [12], A linear relationship 
between model parameters and control voltage is assumed. The MR damper 
hysteresis is reflected in parameters b x and b 2 , while both parameters a v a 0 and b v 
br, model damper viscosity. The force generated by the Bingham damper model, 
F iBD , can be then expressed as follows: 

F w = (, a t + a 2 V )(4 - 4 ) + (&,+ b 2 V) ■ sgn(4 - 4 ) + (c, + c 2 V), (7) 


where: 

c p c 2 -model static parameters, 

V -MR damper input voltage. 

The identification experiment was performed by Plaza in [14], Due to experiment 
results, particularly: symmetrical hysteresis and zero force of piston without load, 
the static property has been simplified, which resulted in Cp c 2 = 0. 

2.3. Inverse MR damper model 

Most semi-active control algorithms require MR damper inverse modelling. An 
inverse model to the Bingham model presented in the previous section can be 
expressed as: 


Vi(F WD , 


) F iBD - a \ (4 ~ 4 ) ~ b \ S S n (4 - Ziu ) 

a 2 (4 - 4 )+ b 2 s s n (4 - 4) 


( 8 ) 
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According to (8), the inverse model is uncontrollable when the denominator 
equals to zero. Considering the MR damper input voltage range, the maximum and 
minimum damper viscosities inputs for inverse model are: c ID mjn = a { -b^ c ID max = 
a j + 5 a 2 + b l + 5 b 2 Inverse model can also be identified by a neuro-fuzzy system [20], 

3. CONTROLLER SYNTHESIS 

Fuzzy controller is a linguistic base-oriented system, which allows for synthesis 
of the control signal on the basis of an expert knowledge. The expert knowledge 
is stored in the form of IF-THEN rules. In this paper the Takagi-Sugeno inference system is 
used, with triangle membership functions defined as {N, Z, P} on the input set [-1,1] 
corresponding to “negative”, “zero”, and “positive” linguistic fuzzy inputs description [20], 
Fuzzy inputs are defined as: x a = (z is - ZJ / g v x 2 = 4 / g 2 , xi 3 = z M / #3’ respectively, where 
g v g 2 , g 3 are constants used for scaling. The fuzzy system output is the damper 
coefficient, and it is given in the form of fuzzy singleton functions {VI, V2, V3, 
V4, V5}, what corresponds to sets evenly spread over the range [V min ,V max ]. The 
control scheme presented in this paper for the 7 DoF model includes independent 
FLC for each wheel. The FLC for a single wheel is presented in Fig. 2. The rule 
base is generated similarly to the skyhook idea, and it is presented in Table 1. A 
benefit of the approach presented in this paper is that no inverse modelling is 
required. However, the rule base is constructed under assumption that the higher 
control voltage the higher force generated by the MR damper. This assumption is 
a kind of simplification due to MR damper hysteresis. The most promising 
advantage of FLC is its high resistance to modelling errors [19]. 


i-th piston velocity 



Figure 2. Fuzzy logic controller for the /-th wheel. 


Table I. 

Fuzzy rule base for FLC 
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3.1. Fuzzy inference 

At first step of fuzzy inference, the degree of membership for each fuzzy input is 
calculated, based on fuzzy sets described above. 

The PROD operator is used for the T-norm. The system fuzzy output is defined 
as a zero-order polynomial, and therefore the fuzzy output V ik is given for each rule, 
as in Table 1. For the system described in such way the FLC output can be defined 
as a weighted sum: 
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Figure 3. Input membership functions for fuzzy controller. 


where w k is a fuzzy weight calculated for each k-th rule: 

w k = H N/pt (x n ) ■ n N/J>k (x n ) ■ n N7pl (x i3 ), (10) 

H N// ,k (x a ) is the degree of membership of input to corresponding fuzzy set {N, Z, P}, 
and li NZP k (x a ) n Nyp k(x a )-dK degrees of membership of inputs x a , x a , respectively [20], 
The fuzzy system output is the MR damper control voltage obtained without a 
need of damper inverse modelling. It is generally hard to heuristically find optimal 
fuzzy parameters, therefore evolution or neuro-fuzzy algorithms are often used 
(e.g. genetic algorithms [7], ANFIS [13]). Those methods can be used offline or 
online. However, their analysis exceeds the scope of this paper. Initial parameters of 
fuzzy controller are based here on the literature, [20], and they can be improved 
according to the user preferences. For example, a stiffer suspension can be achieved 
by narrowing the fuzzy input sets. 

4. SIMULATION RESULTS 

The simulation scenario is presented in Fig. 4. A flat obstacle in the form of a single 
set-off is considered. Vehicle data used for the model come from an all-terrain 

__ fcp 

vehicle presented in [10], They are: M = 670 kg, m i = 30 kg, k iS = 17500 — , 

k iT = 175500-^-, c iD = 1460 — , J FR =1100^-,^ =800% /,fc= 1.9975m, 
nt m m m 

l,r = 0.8025 m Fuzzy logic control simulations are compared th those obtained 
by the Skyhook algorithm in terms of driving comfort and road holding. The sum of 
RMS values of pitch angle and RMS values of roll angle are used for comparison to 
reflect the driving comfort (the higher the value, the worse the driving comfort). The 
comparison of performance of both algorithms reflecting the road holdings 



Figure 4. Simulation scenario (view from the top). 
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evaluated by the RMS value of vertical vehicle displacement with respect to the 
road profile (the higher the value, the worse the road holding). 

The Skyhook damping coefficient in range of c sky = [1000, 10000] is used for 
all comparisons. Obtained results are presented in Figs. 5-7. Scaling parameters 
gj, g 2 , g 3 have been set based on observation of the fuzzy inputs. The choice of 



Figure 5. Vertical displacement of the front-left wheel. 



Figure 6. Vertical displacement of the front-left corner. 




— FLC 

- Skyhook 2000 



Time (s) 

Figure 7. Vertical acceleration of the vehicle centre of gravity (top picture) and vehicle pitch angle (bottom picture). 
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Table 2. 

Comparison between FLC and the Skyhook algorithm (the closer to zero, the better quality indicator) 


Algorithm 

Driving comfort (RMS • I0 3 )[rad] 

Road holding (RMS-10 3 ) [m] 

Skyhook c Sky=|0000 

0.9227 

1.6885 

Skyhook c $ky _ 6000 

0.9327 

1.5927 

Skyhook c Sky = 2000 

1.0657 

1.3300 

Skyhook c S(cy= | 000 

1.2193 

1.2738 

FLC 

1.1279 

1.3774 


appropriate values for these parameters enables to fully utilize the fuzzy linguistic 
rule base. For real applications those parameters may be set based on experimental 
knowledge, vehicle construction analysis, or road profile in-advance observation. 

The quality indicators were calculated for the time between 1 st and 4 th second, 
what corresponds to the time between 1 st impact of the front wheels with the 
obstacle to the last bounce of the rear wheels after 2 nd impact. Namely, the road 
holding indicator stands for the degree of contact between the wheel and the road 
profile. Driving comfort indicator is expressed as oscillations of the vehicle body. 
Both indicators describe the driving safety and passengers comfort. 

The FLC achieved similar road holding and driving comfort indicators, 
comparing to Skyhook with c sky = 2000. However, the FLC does not use any damper 
inverse modelling, in contrast to the Skyhook algorithm. The comparison between 
different c sky values is presented in Table 2. The tuning of FLC has been carried out 
by adjusting scaling parameters. It can also be done by changing fuzzy sets, rule 
base, or fuzzy operators. This is, however, beyond the scope of this paper. 

5. CONCLUSIONS 

This paper shows that fuzzy logic controller used for vehicle suspension vibration 
damping can perform, in some cases, as well as the Skyhook algorithm, but without 
necessity of using an inverse damper model. The use of the fuzzy algorithm is more 
computationally complex than the Skyhook algorithm. However, considering 
cutting edge technology it can be successfully implemented, although with a higher- 
cost controller. There are FLC operations that can be performed once in order to 
reduce the system complexity. Moreover, hardware computing can be used, e.g. the 
Field Programmable Gate Array technology. 

The ease of linguistic description for semi-active vibration control is a promising 
advantage and it deserves more attention. Simulation analysis has been performed 
based on data representing an all-terrain vehicle with four independently suspended 
wheels. 

ACKNOWLEDGEMENTS 

The work reported in this paper has been partially financed by the National Science 
Centre, decision no. DEC-2011/01/B/ST7/06027. 

The authors would like to greatly thank the anonymous reviewers for their 
valuable comments, which helped to improve the paper. 

REFERENCES 

[1] Sammier D, et al.:Skyhook and Control of Semi-Active Suspensions: Some 
Practical Aspects, Vehicle System Dynamics: International Journal of Vehicle 
Mechanics and Mobility, 2003, Vol. 39, No. 4, pp. 279-308. 

[2] Bica I„ et al .'.Physical Characteristics of Magnetorheological Suspensions 
and Their Applications, Journal of Industrial and Engineering Chemistry, 
2013, 19, pp. 394^106. 


Vol. 32 No. 3 2013 





Fuzzy Control for Semi-Active Vehicle Suspension 


[3] Sulaiman S., et al.: Groundhook Control of Semi-Active Suspension For 
Heavy Vehicle, International Journal of Reasearch in Engineering and 
Technology, 2012, Vol. 1, No. 3, pp. 146-152. 

[4] Kaleemullah M., et al.: Design of Robust Hinf Fuzzy and LQR Controller for 
Active Suspension of a Quarter Car Model, Mechatronics (ICOM), 4th 
International Conference On IEEE, 2011, pp. 1-6. 

[5] Tang C., et al.: The Application of Fuzzy Logic Control Algorithm to Active 
Suspensions of Vehicles, Applied Mathematics & Information Sciences, An 
International Journal, 2012, Vol. 6, pp. 855-862. 

[6] Zapateiro M., et al.: Vibration Control of a Class of Semiactive Suspension 
System Using Neural Network and Backstepping Techniques, Mechanical 
Systems and Signal Processing, 2008, Vol. 23, pp. 1946-1953. 

[7] Zhang J., et al.: Fuzzy Logic Controller Based Genetic Algorithm for Semi- 
Active Suspension, Journal of Scientific & Industrial Research, 2012, Vol. 71, 
pp. 521-527. 

[8] Mazur K., Pawelczyk M., Hammerstein Nonlinear Active Noise Control with 
the Filtered-Error LMS Algorithm, Archives of Acoustics, 2013,Vol. 38, 2, 
pp. 197-203. 

[9] Priyandoko G., et al.: Vehicle Active Suspension System Using Skyhook 
Adaptive Neuro Active Force Control, Mechanical Systems and Signal 
Processing, 2008, Vol. 23, pp. 855-868. 

[10] Zareh S. H., et al.: Semi-Active Vibration Control of an Eleven Degrees of 
Freedom Suspension System using Neuro Inverse Model of 
Magnetorheological Dampers, Journal of Mechanical Science and 
Technology, 2012, Vol. 26, No. 8, pp. 2459-2467. 

[11] Kurczyk S., PawelczykM.: Fuzzy Control for Semi-Active Vehicle Suspension, 
Proceedings of the Active Noise and Vibration Control Methods - MARDiH, 
2013, Rytro. - EXTENDED ABSTRACT. 

[12] Krauze P, et al.: Skyhook Control of Front and Rear Magnetorheological 
Vehicle Suspension, In Proc. of 13th Int. PhD Workshop OWD2011, Wisla, 
pp. 380-385. 

[13] Jang J.-S. R., ANFIS: Adaptive-Network-Based Fuzzy Inference System, IEEE 
Transactions on Systems, Man, and Cybernetics, 1993, Vol. 23, No. 3, pp. 665-685. 

[14] Plaza K.: Modelling and Control for Semi-Active Vibration Damping, PhD 
dissertation, Silesian University of Technology, 2008. 

[15] Krauze P: Modelling and Identification of Magnetorheological Vehicle 
Suspension, 17th International Conference on. IEEE Methods and Models in 
Automation and Robotics, 2012, pp. 228-233. 

[16] Kowal J., et al.: Energy Recovering in Active Vibration Isolation System - 
Results of Experimental Research, Journal of Vibration and Control, 2008,14, 
pp. 1075-1088. 


224 


JOURNAL OF LOW FREQUENCY NOISE, VIBRATION AND ACTIVE CONTROL 


Sebastian Kurczyk and Marek Pawefczyk 


[17] Snamina J., et al.: Analiza Rozpraszania Energii Drgariw Aktywnych 
Zawieszeniach Pojazdu dla Wybranych Algorytmow Sterowania, Czasopismo 
Techniczne/Politechnika Krakowska; ISSN 0011-4561. Seria: Mechanika, 
2011, Vol. 108, No. 2, pp. 233-240. 

[18] Sapinski B ., et al.: Experimental Study of a Self-Powered and Sensing MR- 
damper-Based Vibration Control System., Smart Materials and Structers, 2011 

[19] Pawelczyk M.: Integrating Different Control Modes in a Fuzzy Controller, 
Proceedings of 4th IEEE International Conference on Fuzzy Systems, 
Yokohama, Japan, 1995, Vols. I-IV, pp. 1435-1440. 

[20] Cordon O., et al: Applicability of the Fuzzy Operators in the Design of Fuzzy 
Logic Controllers, Fuzzy Sets and Systems, 1997, Vol. 86, No 1, pp. 15—41. 

[21] Pawelczyk M., Kurczyk S.: A Fuzzy Inverse MR Damper Model for Vehicle 
Vibration Control, Proceedings of 20 th International Congress on Sound and 
Vibration, 2013, Bangkok. 


Vol. 32 No. 3 2013 


